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We present a study of the electronic and magnetic properties of the multiple-decker sandwich 
nanowires (CP — M) composed of cyclopentadienyl (CP) rings and 3d transition metal atoms 
(M=Ti to Ni) using first-principles techniques. We demonstrate using Density Functional Theory 
that structural relaxation play an important role in determining the magnetic ground-state of the 
system. Notably, the computed magnetic moment is zero in CP — Mn, while in CP — V a significant 
turn-up in magnetic moment is evidenced. Two compounds show a half-metallic ferromagnetic 
ground state CP — Fe/Cr with a gap within minority/majority spin channel. In order to study the 
effect of electronic correlations upon the half-metallic ground states in CP — Cr, we introduce a 
simplified three-bands Ifubbard model which is solved within the Variational Cluster Approach. We 
discuss the results as a function of size of the reference cluster and the strength of average Coulomb 
U and exchange J parameters. Our results demonstrate that for the range of studied parameters 
U — 2 — AeV and J — 0.6 — 1.2eV the half-metallic character is not maintained in the presence of 
local Coulomb interactions. 
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I. INTRODUCTION 

Utilizing the spin degree of freedom of electrons in the 
solid state electronics has led to the emergence of the 
rapidly developing field of spintronics-. Important elec- 
tronic applications based on the magnetoresitive effect in 
two dimensional heterostructures already exist. Typical 
devices are magnetic read heads and non-volatile mag- 
netic random access memories, where the relative align- 
ment of the layer magnetizations causes large variations 
of the resistance of the structure. This effect is known 
as giant magnetoresistance^, and was discovered in FeCr 
multilayer'^, in which the magnetization of the layers cou- 
ples by the indirect exchange interactions, mediated by 
the electrons of the intermediate layer—. 

There are continuing efforts to improve materials fab- 
rications and design devices for layered magnetic struc- 
tures. With the development of nanotechnology, quan- 
tum structures with dimension of the order of molecu- 
lar or atomic size becomes accessible. In this context, 
the first-principle studies of atomic chain structures, that 
can produce high polarization effects are important. For 
instance, complete spin-polarization in the absence of 
magnetic field was predicted for several bulk materials 
in the class of half-metallic ferromagnets^. The same 
properties were predicted for organic nanowires, such 
as {CeHQ)2V''^ in (C5H^)Fe^ multiple decker nanowires 
or in transition metal doped silicon nanostructures^iiS. 
Organic materials, when used as components of spin- 
tronic devices have significant advantages over inorganic 
ones. First, the spin-orbit and hyperfine interactions 
are weakii, leading to considerably long spin relax- 
ation length and spin-lifetim o ^ ^ i ^ . Organic materials are 



cheap, low-weight, mechanically flexible, and chemically 
inactive, therefore organic spintronics received a consid- 
erable attention, in particular because the spin-polarized 
signal can be mediated and controlled in this case by 
organic molecules^*. 

Metallocene molecules with the general formula 
{Cc,H^)2M consist of two cyclopentadienyl anions [CP] 
bound to a metal center. The metal atom is typically a 
transition metal element from the middle of the 3d-row 
(i.e Ti, V, Cr, Mn, Fe, Co, Ni) as well as from the col- 
umn of iron (i.e. Ru and Os). Molecular structures such 
as CP2V, CP2Cr, CP2Mn, CP2Fe, CP2C0, CPiNi 
has been studied both experimentallyii^"— and theoret- 
ically22. Studies were also reported on the multi-decker 
sandwich clusters based on C^H^ ligands such as triple- 
decker complexes CP^Ni^ and CP^Fe ^^^^^^ and multi- 
ple decker sandwich complexes of type TM{CP2Fe)n+i 
with TM^V, Ti and n = 1 - 32^. The lack of individ- 
ual bonds between the carbon atoms of the CP ring and 
the metal ion makes the CP rings to rotate freely about 
the symmetry axis of the molecule. Although the MT's 
equilibrium structure has a D^d symmetry, the has 
almost the same probability to occur—. Within the equi- 
librium symmetry, d orbitals are splited into a d^^ 
and two sets of doubly degenerate orbitals dxy/dx2_y2 
and dxz/dyz, respectively. 

The electronic structure calculations based on Density 
Functional Theory predict for some metallocene wires a 
half-metallic ferromagnetic ground state — iS^, stimulat- 
ing further the search for such candidate materials for 
spintronic applications^i^i^"— . For all transition metal 
elements the metallocene wires share the same shape de- 
picted in Fig. [TJ It is expected that the p — d covalent 
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FIG. 1: (color online) Geometrical structure of the double 
decker {C5Hz)2M2 sandwich. For the tetragonal unit cell we 
used fixed values for Lx = Ly =12 A while the values of Lz 
are transition metal M-dependent. D illustrates the distance 
between carbon atom and the transition metal atom. 



interactions that take place in metallocene wires are es- 
sential for the physical and chemical properties of these 
compounds. Having similar structure, differences in the 
electronic properties of (6*5^^5)2^/2 can be traced to the 
different electronic filling of the d orbitals that determines 
different magnetic interactions. As d orbitals form rela- 
tively narrow bands, in which electronic correlations are 
essential it is natural to raise the question concerning the 
role of electron-electron interactions in such compounds. 

In the present paper we explore the magnetic prop- 
erties of organometallic multi-decker sandwich wires of 
type {Cc,H^)M where M is a transition metal with atomic 
numbers ranging from Z = 22 (Ti) to Z = 28 (Ni). For 
the interesting case of Cp — Cr in which half-metalicity 
is predicted by Density Functional Theory-based cal- 
culations we formulate a low energy three-band Hub- 
bard Hamiltonian which is solved using many-body tech- 
niques. Our main interest is to understand whether half- 
metallic magnetism is stable against the electronic corre- 
lations in this compound. 



II. COMPUTATIONAL DETAILS 

We have calculated the electronic and structural prop- 
erties of the CP — M wires within the Density Func- 
tional Theory (DFT) using the SIESTA code^^, that em- 
ploys pseudopotentials, and expands the wave functions 
of valence electrons by linear combinations of numerical 
atomic orbitals. We have used a double-zeta polarized 
basis set with an energy cutoff of 30 meV. We have also 



used a number of 16 k-points to model the periodicity 
along the axis of the wire and a 2x2 grid in the trans- 
verse direction. All calculations were performed by us- 
ing generalized gradient approximation to the exchange 
and correlation functional, as parameterized by Perdew, 
Burke and Ernzerhof"^^ . In order to investigate the mag- 
netism in these compounds we consider a unit cell as 
shown in Fig. [l]that would naturally allow to study the 
stability of ferro vs. anti-ferromagnetic states. For all 
systems, we explicitly performed the calculations to look 
for the stable ferromagnetic and antiferromagnetic con- 
figurations allowing structural relaxation. The atomic 
coordinates of all atoms in the unit cell, and the cell's 
length were relaxed up to a maximum force of 0.01 eV/A. 
The values of magnetic moments where computed by in- 
tegrating the spin resolved charge density over "muffin- 
tin" spheres centered on the metal atoms. The radius of 
the sphere was taken 2.25 Bohr , representing about 55 
- 60% of the bond-length M-C for all complexes. All 
the calculations were done using the eclipsed configu- 
ration of the wires (i.e. D^h symmetry). In order to 
study the stability of the wires, we compute the bind- 
ing energy of the metal and CP unit, Ai?, defined as: 
Ai? = Eyjire — 2{Ecp + Em), where the total energy of 
the atom in the unit cell, E^ire, and that of cyclopentane 
unit, Ecp, are computed for the relaxed structures and 
Em is the total energy of the metallic atom. 



III. ELECTRONIC STRUCTURE 

Selected parameters of the relaxed wires are given 
in Table HI We comment in the followings upon the 
atomic-number dependence of the metal-carbon bond- 
lengths. We found that the metal - carbon bond-lengths 
in the investigated CP — M wires display a regular trend. 
Large values are realized for the compounds with extreme 
atomic numbers (i.e. 2.29 A in CP — Ti chain and 2.21 
A in CP — Ni^ respectively) while the minimum occurs 
at the intermediate atomic number (i.e. 2.08A in the 
CP — Mn chain). We note that the same trend occurs 
for the relaxed value of the total length, Lz- The compar- 
ison of the computed carbon - metal bondlength shows 
differences ranging between 0.5 % (for CP — Ni ) and 
about 3 % (for CP — V). An important difference occurs 
for CP - Mn (i.e. 15%). We note that for CP - Mn 
we also found a ferromagnetic state which is less sta- 
ble than the non-magnetic one (energy difference is 0.15 
eV); for this state the C-Mn bondlength is 2.12 A. This 
magnitude of the bond-length is in agreement with the 
experimentally reported value determined for a high spin 
state (see Ref. Il5l ). The vanishing of magnetic moment 
in the periodic structure is probably the consequence of 
the large differences between the geometric structure of 
the molecule and the corresponding chain. 

As one can see in Tab. Uthe values for the binding en- 
ergy AE, depend on the gradual filling of the d-orbitals 
of the metal atom as was observed previously^*'. The 
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Atom 


Lz 


D 


Afo 




Type 




[A] 


[A] 






[eV] 



Ni(3d**4s^) 7.35 2.21 (2.20") 0.00 0/0 NM -6.91 
Co(3d^4s^) 7.10 2.17 (2.12*') 0.00 1.31/-1.31 AFM -11.38 
Fe(3d®4s^) 6.77 2.10 (2.06'') 2.00 1.00/1.00 HMF -11.06 



0.00 



0/0 



NM 



Ma{2,dP'is^) 6.59 2.08 (2.38 , 

Cr(3d'^4s^) 6.90 2.14 (2.17") 2.00 1.00/1.00 HMF 

V(3d^4s^) 7.20 2.21 (2.28^) 3.99 1.70/1.70 FM 

Ti(3d^4s2) 7.52 2.29 (2.36 



8.38 
-7.98 
13.73 

0.00 1.36/-1.36 AFM -14.98 



TABLE I: Geometric and magnetic properties of the wires. 
Lz is the length of the unit cell, D is the metal - carbon 
distance (see also Fig. [1]). Mo is the magnetic moment 
per cell, M1/2 are the magnetic moments for the transition 
metal atoms. AFM/FM - antiferro/ferro-magnetic; NM - 
non-magnetic, HMF - halfmetallic ferromagnetic AE is the 
binding energy of the chain. Experimental data on isolated 
molecules are indicated for the metal-carbon bondlength. 
(" Reference2i, Referencei^, " Reference^^, Reference,— 
and^, " and ^ Referencei^, » Reference^i) 



strongest binding energy is realized for CP — Ti; fur- 
ther filling the d-orbital reduces the strength of binding 
(i.e. non-bonding states are populatec^). An almost 
twice weaker binding energy is obtained for CP — Cr (- 
7.98eV) in comparison with Cp — Ti(-14.98eV). Increas- 
ing the number of d electrons starting from Cr Mn 
Fe Co the binding energy is increasing again, while at 
CP — Co, the partial filling of the anti-bonding orbitals 
stops this trend. A further increase of the electronic pop- 
ulation in the anti-bonding orbitals of nanowires, lead to 
an important reduction of the binding energy as seen for 
CP - Ni wire. 

Total densities of states for the seven wires are given in 
Fig. [2I Increasing the number of valence electrons from 
Ti {Z = 22) to Fe {Z = 26) resuhs in a gradual fiUing 
of the d band, therefore the Fermi level is expected to 
be shifted towards higher energies, in addition to which 
the magnetic interactions contribute in determining the 
magnetic ground state. In CP — Ti and CP — V we found 
that both ferromagnetic and antiferromagnetic states are 
possible. For Ti the ground-state AF configuration is 
lower in energy by 0.34eV, while for V the ferromagnetic 
is most stable with an energy difference of about 0.1 7eV. 
CP — Mn shows long-range ferromagnetic order for a 
non-relaxed structure that disappears in favor of a non- 
magnetic ground state in the relaxed case. For all the 
other compounds a single stable solution has been ob- 
tained, for the considered range of parameters and con- 
figurations. Two systems where found to be half-metallic: 
CP - Cr and CP - Fe. For CP - Cr the gap is situ- 
ated in the majority spin DOS and it has a value of 0.48 
eV, while for the case of CP — Fe the gap occurs in the 
minority spin DOS with a magnitude of about 0.41 eV. 

The DOS projected onto the carbon atoms (not shown) 
has similar shapes as the total density of states but with 
a significantly smaller weight. In order to determine the 
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FIG. 2: (color online) Total DOS for the computed wires. 
With red line we show the cases when half-metalicity occurs 
for the CP — Fe (minority spin gap) and CP — Cr (majority 
spin gap) sandwiches. 



character of the orbitals near Fermi level for the interest- 
ing case of half- metallic CP — Fe and CP — Cr we com- 
puted the partial densities of states. In the upper panel of 
Fig. [21 the results for Cp—Fe are seen. A significant con- 
tribution within the majority spin channel are provided 
by the degenerate Fe — 3dxz/-idyz partially occupied or- 
bitals. The same orbitals are completely empty in the 
minority spin channel. The other 'id^-2,'idxy/'idx^_y-2 or- 
bitals are completely occupied, so that the gap is realized 
between the minority spin ^dxy/Sd^^-y^ and Sd^zf-^dyz 
orbitals. For the case of the majority spin half- metallic 
Cp — Cr, the degenerate Cr — Sdxz/Sdyz, are empty in 
both majority and minority spin channels. Partial oc- 
cupation is obtained for minority spin Sdxy / 'id^^ -y^ and 
3^22 orbitals, while the same orbitals are completely filled 
within the majority spin channel. Contrary to Cp — Fe 
in Cp — Cr the gap is realized between the majority 
spin 3dxy/-idx2-y2 and ?>dxz/^dyz orbitals. For both half- 
metals, the carbon's Pz orbitals are also present at the 
Fermi level, however presenting only a weak spectral 
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weight. Note that, DOS projected over the Pz orbitals 
of a carbon atom in CP ring was multiphed by five to fa- 
cihtate the direct comparison with the total DOS shown 
in Figl21 The computed magnetic moment per cell for 
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FIG. 3: (color online) Partial density of states for the half- 
metallic CP — Fe and CP — Cr wires. The degenerate M-Scia;^ 
and yi-Zdyz orbitals are represented with blue lines, red line 
represents yi-Zd^z and }A-<id^2_y2, and with black line the M- 
'id^2 orbitals are shown (M=Fe,Cr). Finally with green line 
we show the sum over the five C-pz orbitals. 

the investigated wires takes values between and 3.99 
^B, and can be attributed mainly to the transition met- 
als, while the CP rings have very weak magnetic mo- 
ments. The expected monotonic increase of the magnetic 
moment as function of d-orbitals occupancy can be ex- 
plained by considering that the metal atom is donating 
one electron to the CP ring, which forms a six electrons 
TT ring. Therefore, N^ai — 1 electrons rest on the metal 
and are determinant for its magnetic state, where N^ai 
is the number of valence electrons. By using Hund's rule 
and the transition metal valence electronic configuration 
(i.e. 4s, 3^22, 3dxy/3dx2-y2 and Sdyz/Sdzx: respectively) 
we can assign a magnetic moment to each wire. We ex- 
pect therefore the magnetic moments of 0,2,4,2,0 and 2 
for CP - Co, CP - Fe, CP - Mn, CP - Cr, CP~V 
and CP — Ti, respectively. Indeed, the electronic struc- 
ture calculation provides these results with two remark- 
able exceptions: CP — Mn {0 fiB - instead of 4 (Xb) 
and CP — V (3.99 - instead of /z_b)- Note also 
that Cp — Ti is an exception from the above discussed 
Hund's rule as its ground state is anti-ferromagnetic in- 
sulator with moments of ±1.36. According to our results. 



the structural relaxation effects determines the departure 
from the expected behavior in all three compounds: the 
significant drop/turn- up of the magnetic moment in the 
CP — Mn/V and the weaker reduction of the magnetic 
moment in the case of Cp — Ti . We have additionally 
checked that Mn replacing Cr in the frozen CP ~ Cr 
structure, retains a non-zero magnetic moment (i.e. 2.36 
jiB ), with values close to the Cr ones (i.e. 2.0 ^b)- An 
overall integer magnetic moment on the unit cell (the 
value Mq shown in table Table H]) is obtained for all the 
compounds. 




FIG. 4: (color online) Real-space spin density contour plots 
for CP — Fe. Majority/minority spin density is represented 
with blue/red lines. The black lines represent the regions 
with no difference between up and down spin densities. Up- 
per/lower picture presents a side/top view (see text for geo- 
metrical details of the contours). 

In Fig|3] we illustrate with two contour plots the dif- 
ference between majority and minority real space spin 
densities, n(r)''' — n(r)^ for the Cp-Fe wire. We use two 
rectangular planes to plot the contours. The first plane 
has the edges parallel to the z axis and is passing through 
the metal atoms and one of the carbon atoms in the CP- 
ring (top panel Fig. H]). This plane provides us a side 
view, with respect to the axis of the nanowire. The lower 
panel of Fig. U illustrates the "second" plane that has 
the edges parallel to the x and y axes and provides a 
top view of the spin density in the cyclopentane ring. 
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For both spin orientations we use ten equidistant con- 
tours and an increment of ± 0.0025 e/Bohr^ for the up 
/ down spin densities. Note that the spin density on car- 
bon atom has a maximum value close to 0.0075 e/ Bohr^ 
(i.e. three contour lines). As seen on the contour plot, 
this maximum value is at least two orders of magnitude 
lower than the values obtained for the spin densities on 
the metal site. For Fe the maximum spin density is lo- 
cated around its atomic position and has a value of about 
1.2 e/Bohr^ (picture not shown here). In addition for the 
Fe sites the difference of real-space spin densities has a 
positive sign (blue contours), while on the carbon sites 
the difference is negative signaling an anti-parallel orien- 
tation of a small magnetic moment with respect to Fe 
magnetic moment. 

For CP — V a, non-zero spin-density contribution is 
present on the CP-rings, which makes this material to be 
an almost half metallic compound. We also noticed that 
in the case of CP — Co, the majority spins accumulate on 
one side of the cyclopentane ring, therefore contributing 
to the anti-ferromagnetic properties obtained for CP — 
Co. 



IV. MULTI-ORBITAL HUBBARD MODEL FOR 
HALF-METALLIC CP-CR CHAIN 



symmetrically orthonormalized NMTOs is a set of Wan- 
nier functions. In the construction of the NMTO basis 
set the active channels are forced to be localized onto 
the eigenchannel RZm, therefore the NMTO basis set 
is strongly localized. For Cp — Cr we choose to down- 
fold (i.e. to integrate out) all orbitals except the Cr-3d 
{xy,x^ —y^jZ^ — l) manifold. At this point is important 
to mention that the NMTO calculations has been per- 
formed considering the structure shown in Fig. [T] with 
the parameters discussed in Table HI The calculations 
considered the tetragonal symmetry in such a case the 
Cr-3d(x2/,x^ — y^) orbitals are nearly degenerate, while 
in the computation using SIESTA a complete degener- 
acy is seen as the D^h structural symmetry of the ligand 
field is considered. This difference turns out not to be 
essential in considering the many-body effects presented 
here. Thus, the effective Hamiltonian is confined to a set 
of effective three orbitals of {xy, — y^, — l)-type in 
a reduced window of energies. For optimizing the energy 
window with respect to the orbitals we chose the fol- 
lowing expansion points with respect to the Fermi level: 
Ey - Ef ^ 0.31el/, -1.05eF and -1.59eT^, where the 
Fermi energy has the absolute value Ep = — 3.44eV. In 
Fig. [5] we show the eigenvalues of the effective Hamilto- 
nian along some high-symmetry directions in comparison 
with the full orbital basis. 



A. Model setup 

In order to discuss further the magnetic properties 
of the half-metallic Cp — Cr metallocene nano-wire we 
supplement the DFT analysis with many-body calcula- 
tions based on a multi-orbital Hubbard model. Such 
calculations were performed recently for the prototype 
half-metallic ferromagnet NiMnSb including local'^^ and 
non-local correlation'^^''^^ effects, demonstrating the ap- 
pearance of significant many-body induced states within 
the half-metallic gap leading to the disappearance of 
half-metalicity. Our primary goal is to investigate the 
existence of similar effects in the half-metallic met- 
alocene nano- wires, in particular the stability of the half- 
metalicity in the presence of the Coulomb interaction. 
In the present section we discuss the construction of the 
low energy Hamiltonian in a similar way as has been re- 



cently performed for NiMnSb 



36,37 



CrO 



38 



or TiN^ 



the framework of N-th order muffin-tin-orbital method. 
Alternative downfolding schemes, can be also formulated 
within tight-binding approaches^Siii. 

The starting point in the combined electronic struc- 
ture and many-body calculation is the formulation of the 
low energy model Hamiltonian. Specifically, the uncorre- 
lated part of the Hamiltonian is obtained from the down- 
folding procedure^i^ within the N-th order muffin-tin- 
orbital (NMTO) method. The NMTO raeiho^B^ can 
be used to generate truly minimal basis sets with a mas- 
sive downfolding technique. Downfolding produces bands 
with a reduced basis which follow exactly the bands ob- 
tained with the full basis set. The truly minimal set of 
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FIG. 5: (color online) Band-structure of CP — Cr around 
Fermi-energy obtained with the LDA calculation (red line) 
for the complete orbital basis set and the NMTO bands ob- 
tained after downfolding to the effective xy, —y^ and 2^ — 1 
orbitals. 

Fourier-transformation of the orthonormalized NMTO 
Hamiltonian, iJ'"^^(k), yields on-site energies and hop- 
ping integrals, 



ttLDA _ / _L \ajLDA _ I ± \ _ .a: 
-"Om'.Rm — \AOm' | \ ARm/ — ''n 



_ ,xyz 



(1) 



in a Wannier representation, with the orthonormal 
|XR,m)-NMTO Wannier functions, m is labeled by the 
three xy, x^ — y^, — \ orbitals. In this representa- 
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tion the on-site matrix elements „ are nearly diago- 
nal, a very small non-zero coupling between the — 
and — \ orbitals is obtained, which is a consequence 
of the considered tetragonal symmetry. A set of rotated 
orbitals has to be introduced such that the local single- 
particle density matrix is diagonalized. Such a procedure 
is frequently used'*^ and in the present case this leads to 
a small change in the effective hopping parameters of 
Cp — Cr. The on-site and the directional hopping matrix 
elements within the rotated basis-set are given (in units 
of eV) by: 



.000 



.00±1 



^ -4.245 

-4.287 
-3.869 

^ -0.545 

-0.589 -0.02 
-0.02 0.091 



(2) 



Thus, the non-interacting part of the effective Hamilto- 
nian for CP — Cr has the form: 

E ^m'r™ 4™vcr™. + /i-c- (3) 

R,' ,R,,{m' ,m},cj 

To take into account correlation effects to the non- 
interacting Hamiltonian in eq. [3] we add the Hubbard 
part such that a 3-band "correlated" Hamiltonian is cast 
in the form: 



R,m 

R,,m<m' ,(7,cj' 
I{..7n<m' 

+ ^ '^mm'cl^^^f^cl^^'^C^rnicl^^^-^ ^ h.C. (4) 
R,,m<m' 

where Hq is given by eq. (O, nR„„ = clima^Rma is 
the number of particle operator and c^mai'^^""^) 
the usual fermionic creation(annihilation) operators act- 
ing on an electron with spin a at a site R in the or- 
bital TO. The on-site Coulomb interactions are expressed 
in terms of two parameters U and J via: Umm = U, 
Umm'i^m) ^ U' ^ U ~ 2J, J,„„/ = Jii. In our cal- 
culations we used values of U in the range of 2 — AeV 
and J in the range of 0.6 — 1.2eV as they are usual pa- 
rameters for transition metal elements^^— . Similar val- 
ues were used to include correlation effects within mean- 
field GGA-I-U calculations^Si^ in multiple-decker-type of 
compounds. Note that corrections for a double-counting 
of the interaction within the manifold of orbitals consid- 
ered, only produces an irrelevant constant shift of the 
chemical potentiaP^'^l as we are considering a model 
Hamiltonian with fixed number of electrons. 



B. The correlated ferromagnetic state 

The total Hamiltonian in eq. is solved by the VGA 
method, that is an extension of Gluster Perturbation 
Theory (GPT)S.-^. In this approach, the original lat- 
tice is divided into a set of disconnected clusters and the 
inter-cluster hopping terms are treated perturbatively. 
VGA additionally includes "virtual" single-particle terms 
to the cluster Hamiltonian, yielding a so-called reference 
system, and then subtracts these terms perturbatively. 
The "optimal" value for these variational parameters is 
determined in the framework of the Self-energy Func- 
tional Approach (SFA)^^'^^, by requiring that the SPA 
grand-canonical potential f2 is stationary within this set 
of variational parameters. 

- CP-Cr-CP-Cr-CP-Cr-CP-Cr 




Reference system 
2-sites cluster 



Reference system | Reference system 
3-sites cluster 4-sites cluster 



FIG. 6: (color online) Schematic representation for the model 
CP — M chain. Three different choices for the reference sys- 
tems are shown: the 2,3 and the 4 sites clusters. Intra- and 
inter-site couplings are described by the t^~^ matrix, with e 
the on-site energies for the effective d-orbitals rotated in the 
local basis. 



In fig. [S]we show the choice of the reference system, ob- 
tained by disconnecting the multiple decker nanowire (i.e. 
the chain of CP — M units) into a set of sites connected 
by the same intra- and inter-cluster hopping t^,~^- The 
electronic states connected to the CP-rings and some of 
the lower and higher d— energies were integrated-out so 
that the electronic structure around the Permi level is 
described by the low energy Hamiltonian Eq. 2] formed 
with the three active effective orbitals. 

Our finite size scaling analysis has been performed for 
the fixed values oi U — 2eV and J = 0.6el^ parameters 
for the different clusters of two, three and four sites. The 
results indicate that larger the cluster size convergence 
and optimization of the VGA grand potential is still eas- 
ily achieved. In particular, the self-consistent solutions 
for the grand potential, in the case of three and four 
sites clusters provides similar occupations to the non- 
interacting case. The non-interacting results were ob- 
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tained by self-consistent calculations of the grand poten- 
tial in which the exchange splitting has been considered 
as variational parameter parameter. For the four sites 
cluster, a total magnetic moment of 1.10/is is obtained 
within the non-interacting calculation while for the in- 
teracting case a value of about 0.95 — 0.97/iB is obtained 
depending on the strength of average Coulomb parame- 
ters {U and ,/). 



ble peak structure just below Fermi level. Increasing the 
size of the cluster (see left- and right-column of Figl?]), 
there is no significant change within the minority spin 
channel. All essential features of density of states remain 
unchanged when comparing the results of calculations for 
three and four sites clusters. A close look reveals only the 
scale difference connected of the number of sites consid- 
ered within the cluster. 
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FIG. 7: (color online) Orbital resolved densities of states of 
CP — Cr chain for the clusters containing three (left column) 
and four (right column) sites, and the average Coulomb and 
exchange interactions are U — 2eV and J — O.GeV. At the 
non-interacting level (dashed-blue line) a half-metallic solu- 
tion is obtained, while in the presence of interaction (solid- 
red line) majority spin gap is closed and a metallic solution 
is obtained. 

In fig. [7]we show the total and orbital resolved density 
of states, for the clusters of three and four sites and fixed 
values U — 2eV, J = 0.6eV. As one can see, the majority 
spin gap is filled and normal ferromagnetic ground state 
is obtained with a considerable spectral weight just below 
the Fermi level. In comparing to the non-correlated cal- 
culations (dashed-blue line), most significant changes are 
seen in the majority spin channel for all three orbitals. 
The majority — 1 orbital is pushed towards the Fermi 
level and splitted by electronic correlations, nevertheless 
it remains completely occupied. The xy and x'^ — 
orbitals remain almost degenerate and are shifted also 
towards the Fermi level. In addition they develop a dou- 
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FIG. 8: (color online) Orbital resolved densities of states 
for the four site cluster (interacting/non-interacting - red 
solid/blue dashed line) with increasing values of U and J as 
indicated in the figure. 

In Fig. [5] we present the results for the four site clus- 
ter, for two distinct values of U and J. Increasing the 
values of U and J, the density of states features situated 
between —4 and —3eV are further renormalized towards 
higher binding energies. Within an energy window of 
ztleV around Fermi level (see Fig. |9]) there is a slight 
shift of minority spin states towards Fermi level, while 
for the majority spins, spectral weight transfer towards 
Fermi level is slightly increased for larger U. For the ma- 
jority electrons, DOS in the close vicinity of Fermi level 
obtained within the many body calculations shows the 
presence of some states pinned almost at the Fermi level. 
The position of these pinned states remains practically 
unchanged, for the different U and J parameters, while 
their spectral weight is slightly increasing with U and 
J. From the orbital projected densities of states seen 
in Fig. [7] and Fig. [5] we can identify its composition 
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as being mainly determined by the xy and x'^ — or- 
bitals. Around 0.25eV below Fermi level, the peak hav- 
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FIG. 9: (color online) Total density of states for the four 
site cluster for different values of U and J parameters. Inset 
presents the integrated spectral weight of the peak centered 
at Ef- 

ing mainly a — 1 character is seen to shift towards Ep 
when U and J increase. In addition spectral weight is 
transferred towards the peak pinned at the Fermi level. 
Within the minority spin channel, increasing U/J leads 
to a visible shift of the — 1 orbital towards and above 
the Fermi level. The inset of Fig. IH] shows the integrated 
spectral weight of the states pinned in the close vicinity 
of Fermi level. A similar effect has been recently dis- 
cussed for the ferrimagnet Mn2VAl in which the local 
but dynamic electronic correlations captured by DMFT^^ 
demonstrated the closure of the gap and a very strong 
depolarization effect. For both these compounds LDA 
based DFT calculations predict a half-metallic majority 
spin gap, but for both materials correlations effects are 
shown to destroy half-metalicity, by creating many-body 
induced states just below the Fermi level. Because of 
their position below Fermi level and their considerable 
spectral weight, one can expect that these states might be 
detectable by spin-polarized photoemission, which con- 
stitutes an interesting possibility to investigate many- 
body effects in such organo-metallic spintronic materials. 

V. CONCLUSION 

In this work we studied the electronic structure and 
magnetic properties of a series of metallocene multiple- 



decker sandwich nanowires with the formula {CzH^)2M2: 
where M=Ti to Ni. Our results show that a broad vari- 
ety of electronic and magnetic properties are predicted in 
these wires (i.e. insulator, metallic, half- metallic). Based 
on the results of the binding energies, we show that all 
the structures are energetically stable. We demonstrate 
that physical properties of the multiple-decker nanowires 
are determined by structural relaxation. Depending on 
the transition-metal element we found regular alterna- 
tions of the geometric parameters such as bondlengths 
and cell length, that determines the magnetic proper- 
ties of {C^Hc,)2M2 systems. Accordingly we found that 
CP — Co is an anti-ferromagnetic insulator with a large 
gap of about 1.5eV. The Mn moment's reduction to 
zero in CP — Mn and the considerable turn-up of the 
V moment in CP — V is obtained. For the CP — V and 
CP — Ti both ferro and antiferro solutions are possible 
which allows us to estimate the exchange couplings be- 
tween the V-V and Ti-Ti atoms. These are of magnitude 
Jv = OMeV « 3945K and Jn = O.lTeV « 1973ii: and 
the corresponding ground state in these cases is metallic 
ferro and anti-ferromagnetic respectively. In all calcula- 
tions the CP— rings carry a very small induced magnetic 
moment, that contributes into the stabilization of the 
structure and the transmission of the magnetic interac- 
tions through the wire. For the Fe-based and Cr-based 
sandwiches we found within the DFT - GGA calculations 
stable half-metallic solutions. 

For the CP—Cr half- metallic wire we have constructed 
a low-energy Hamiltonian by downfolding all orbitals ex- 
cept Cr-3d{xy, x^ — y^ , — 1) orbitals which provides 
an effective low-energy model. This simplified model- 
Hamiltonian is used to investigate correlation effects in 
the half-metallic ground state. We have solved the three- 
band Hubbard Hamiltonian within the Variational Clus- 
ter Approach. The main effects determined by electronic 
interactions are: (i) within the minority channel (spin 
down) the — 1 orbital is less occupied and more itiner- 
ant, while the other minority spin orbitals do not suffer 
significant change; (ii) for the majority electrons (spin 
up channel) many body state of xy and x^ — origin 
are pinned within the half-metallic gap just below Fermi 
level. Their position is practically unchanged for different 
values of U and J, in contrast to the occupied z^ — 1 or- 
bitals that shifts towards the Fermi level, with increasing 
U . The overall conclusion is that the majority spin half- 
metallic ground state is not robust against the presence 
of local Coulomb correlations, and a correlated ferromag- 
netic ground state is obtained in which the majority spin 
gap is closed by many-body induced states. 
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